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QCD thermodynamics is studied using a model which combines Polyakov loop dynam- 
ics with spontaneous chiral symmetry breaking and its restoration (the PNJL model). The 
input is fixed entirely by pure-gauge lattice QCD results and by pion properties in vac- 
uum. Successful comparisons with results from thermal lattice QCD are achieved, including 
extrapolations to finite quark chemical potential. The phase diagram and selected suscepti- 
bilties for two quark flavors (Nf = 2) are investigated with inclusion of diquark degrees of 
freedom. 

§1. Introduction 

This presentation addresses the following principal question: is it possible to un- 
derstand the results of lattice QCD (LQCD) thermodynamics in terms of properly 
selected quasiparticle degrees of freedom? A suitable framework for approaching this 
question is the PNJL model, a synthesis of Polyakov loop dynamics and the Nambu & 
Jona-Lasinio (NJL) model. The Polyakov loop (or thermal Wilson line) is an indica- 
tor 1 - 1 of the confinement — > deconfinement transition in QCD. The NJL model, 2 - 1 ' 3 ) 
on the other hand, provides a schematic description of the dynamical mechanism 
which drives spontaneous breaking of chiral symmetry in QCD. The PNJL model 
is thus based on two seemingly disconnected symmetries which characterise QCD in 
opposite limiting situations: 

• The Z(3) center symmetry of the SU(3) color gauge group is exact in the limit of 
pure gauge QCD, realised for infinitely heavy quarks. In the high-temperature, 
deconfinement phase of QCD this Z(3) symmetry is spontaneously broken, with 
the Polyakov loop acting as the order parameter. 

• Chiral SU(Nf)ji x SU(Nf)i symmetry is an exact symmetry of QCD with Nf 
massless quark flavors. In the low-temperature (hadronic) phase this symmetry 
is spontaneously broken down to the flavor group SU(Nf)v (the isospin group 
for Nf = 2 and the "eightfold way" for Nf = 3). As consequence there exist 
2Nf + 1 pseudoscalar Nambu-Goldstone bosons and the QCD vacuum is non- 
trivial. It hosts quark condensates (qq) which act as chiral order parameters. 
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Confinement implies spontaneous chiral symmetry breaking, whereas the re- 
verse is not necessarily true. Whether and under which conditions the chiral and 
deconfinement transitions coincide, as apparent in some lattice QCD results, is thus 
a fundamental issue. Numerical simulations of QCD thermodynamics on the lat- 
tice 4 )' 5 ) are the primary sources of information for our purpose. The equation of 
state of strongly interacting matter is now at hand as a function of temperature T 
and in a limited range of quark chemical potentials fi. Strategies for circumventing 
the fermion sign problem characteristic of finite \x are: improved multi-parameter 
re-weighting techniques, 6 ) Taylor series expansion methods 7 ) -9 ) and analytic contin- 
uation from imaginary chemical potential. 4 )' 10 )' n ) LQCD data exist for the pressure, 
the energy and entropy densities, quark densities and selected susceptibilities. 

The PNJL model is a useful device in order to interpret such LQCD results, to 
understand the underlying systematics in terms of quasiparticles and to extrapolate 
into regions not accessible by lattice computations. This model can also contribute to 
the question about the entanglement between chiral and deconfinement transitions. 
It is designed much in analogy with a Ginsburg-Landau type approach: identifying 
the relevant order parameters as collective degrees of freedom which govern the 
dynamics and thermodynamics. Earlier versions of such a model have been discussed 
in refs. 12 )' 13 ) The present report summarises our recent results and developments of 
the two-flavor PNJL model 14 ) -17 ) including successful comparisons with a variety of 
LQCD data. For related works see refs. 18 ) -20 ) 



§2. Sketch of the PNJL Model 

The two-flavor PNJL model 16 ) is specified by the Euclidean action 



P=l/T 



dr 



d 3 x 



v 



^d T ^ + H(ip,^,(f))\+dS{<f),T) (2-1) 



with the fermionic Hamiltonian density*); 

H = —iip^ (a • V + 74 mo — 4>) tp + V(tp, ifi) 



(2-2) 



where ip is the Nf = 2 doublet quark field and mo = diag(m u ,md) is the quark 
mass matrix. The interaction V(ip,^), to be specified in detail later, includes four- 
fermion couplings acting in quark- antiquark and diquark channels. The quarks move 
in a background color gauge field (p = = iAq , where Aq = gAa t a with the 
SU(3) C gauge fields A% and the generators t a = X a /2. The matrix valued, constant 
field <p relates to the (traced) Polyakov loop as follows: 



V exp 



i I drA^ 



-Tre** 
3 



(2-3) 



In a convenient gauge one can choose a diagonal representation, 4> = ^3 ^3 + 08 ^8 , 
which leaves only c/^g as field variables representing the Polyakov loop. 



*) 



a — 70 7 and 74 = 470 in terms of the standard Dirac 7 matrices. 
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2.1. Polyakov loop effective potential 

The piece SS = -^U((p,T) in (2-1) carries information about the gluon dy- 
namics which, in the present approach, is approximated by the dynamics of the 
Polyakov loop. The effective potential IA models the region around the confinement- 
deconfinement transition in pure gauge QCD on the mean field, Ginsburg-Landau 
level. One can expect such an approach to work up to temperatures around twice 
T c . At much higher temperatures a description of the thermodynamics entirely in 
terms of the Polyakov loop is no longer adequate as transverse gluons will become 
important. 

The construction of the Landau effective potential describing Polyakov loop 
dynamics is guided by the Z(3) center symmetry which transforms an element 
u £ SU(3) C to exp(27rm/3)it (n = 1,2,3,...). The basic building blocks for such 
a potential are <P*<P, <P 3 and <P* ? ' terms. In an earlier version of the model 14 ) a poly- 
nomial form for IA was introduced and subsequently used also by other authors. 18 ) 
However, this form does not satisfy the correct high-temperature limit. In the up- 
dated version 15 ) -17 ) the effective potential IA has been designed to properly meet 
group theoretical constraints, with an ansatz motivated by the SU(3) Haar measure: 



(2-4) 



£/(<l>,^> , T) = -l a (T) <p*$ + b(T) In 1 - 6<2>*<2> + 4 (V 3 + <2> 3 ) - 3 {$*$f 

The temperature dependent prefactors are written 

a(T) =a + ai + a 2 ' , b(T) = b 3 ' . (2-5) 

The logarithmic divergence in Eq.(2-4) near <P*, <P — > 1 makes sure that the Polyakov 
loop does not exceed constraints for the normalized trace of an SU(3) element. The 
parameters are chosen such that the pressure p = —U (T) and related thermodynam- 
ical quantities reproduce the LQCD results for the pure glue equation of state 25 ) 
with its first order phase transition at a critical temperature To ~ 270 MeV. The 
numerical values are given in ref. 15 )' 16 ) : 

a = 3.51 , ai = -2.47 , a 2 = 15.2 , b 3 = -1.75 . 

Uncertainties in ai 5 2,3 are estimated to be at the level of a few %, while ao = 
by virtue of the Stefan-Boltzmann limit. Figs. 1,2 show the resulting fit to the pure 
gauge QCD equation of state and the Polyakov loop effective potential displaying 
the first order transition at T = Tq. 

2.2. NJL effective interaction 

The interaction V in Eq.(2-2) includes chiral SU(2) x SU(2) invariant four-point 
couplings of the quarks acting in quark-antiquark and diquark channels: 



(^Y + (t/u7 5 tV) 2 - ^ pC 75 T 2 A 2 ^ T ) (V> T 75T2A 2 C^)] , (2-6) 



H 
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Fig. 1. Scaled energy density e, entropy density 
s and pressure p as functions of temperature 
(given in units of the critical temperature 
T = 270 MeV) of pure gauge QCD. Lattice 
results 21 ' are compared with the Polyakov 
loop model (curves) 16 ' with parameters as 
explained in the text. 



Fig. 2. Effective potential 2-4 for the Polyakov 
loop $ at different temperatures (given in 
units of the critical temperature To = 270 
MeV) with parameters adjusted to repro- 
duce the pure glue QCD equation of state 
(see Fig.l). 



where C is the charge conjugation operator. One can think of Eq.(2-6) as a subset 
in the chain of terms generated by Fierz-transforming a local color current-current 
interaction between quarks, 

In this case the coupling strengths in the quark-antiquark and diquark sectors are 
related by G = ^H, the choice we adopt. The minimal ansatz (2-6) for V is motivated 
by the fact that spontaneous chiral symmetry breaking is driven by the first term 
while the second term induces diquark condensation at sufficiently large chemical 
potential of the quarks. Additional pieces representing vector and axialvector qq 
excitations as well as color-octet diquark and qq modes are omitted here. Their 
effects are not important in the present context. 

The NJL-model with Nf = 2 involves three parameters: the current quark mass 
m u d, a local four -quark coupling strength G and a three-momentum cutoff A. As 
in refs. 14 )~ 16 ) we choose 

m u ,d = m = 5.5 MeV , G = -H = 10.1 GeV -2 , /L = 0.65 GeV. 

3 

These values were fixed to reproduce the pion mass and decay constant in vacuum, 
m n = 139.3 MeV, f n = 92.3 MeV, and the chiral condensate, {ipu^u) = —(251 
MeV) 3 . 

§3. Thermodynamics of the PNJL Model 

In order to evaluate the thermodynamical properties of the PNJL model, new 
auxiliary fields (a, tt) and [A, A*) are introduced by bosonisation, absorbing the 
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relevant quark-antiquark and diquark (antidiquark) correlations, respectively. The 
bosonised action involves the trace-log of the inverse Nambu-Gorkov matrix propa- 
gator which incorporates the couplings of quark quasiparticles, scalar and diquark 
fields. The pion field does not acquire an expectation value and the thermodynamic 
potential at mean field level becomes, after carrying out Matsubara frequency sums: 



f2(T,n)=U($, T) + — 4 

-2N, 



a 2 A*A 



2G 2H 

d 3 p 



(2*)' 



1 + e" 



-Ej/T 



(3-1) 



The quasiparticle energies Ej are determined analytically, with the results: 



Ei, 2 = e(p)T (m + 



2i> 



V3 



3,. ..,6 



e{p) ± 



II 



V3 



+ |Z\| 2 =F i ' 



>3 ■ 



Here e(p) = \f p 2 + m 2 . The mass of the quark quasiparticles is generated dynam- 
ically by the gap equation, m = niQ 4 a = tjiq — Gtyip). The energy difference 
AEj = Ej — £o i n is taken with respect the free fermion energy, eq = \Jp 2 + mg. 

At non-zero chemical potential fj,, the effective action is complex. The fermion 
sign problem persists in the PNJL model just as in QCD itself. As shown in ref. 22 ) 
and following discussions in ref., 23 ) the consistent minimisation condition for the 
thermodynamic potential is 



dRefl 



d (a, A, fa, fa) 



. 



(3-2) 



In general, the Polyakov loop fields <P and <P* (and their thermal expectation val- 
ues) are different at non-zero quark chemical potential. This is primarily an effect 
of fluctuations of the fields around their mean values. The mean field approxima- 
tion identifies the fields with their expectation values and implies 22 ) <P = or 
equivalently, ^8=0. 



§4. Results 



We now proceed to discuss selected examples and applications of PNJL thermo- 
dynamics, mostly in comparison with thermal LQCD results. 

4.1. Chiral condensate and Polyakov loop 

Consider first the chiral condensate {tpip) and the Polyakov loop <P and their 
temperature dependence at zero quark chemical potential (see Fig. 3). The original 
first-order deconfinement transition of pure gauge QCD is now turned into a crossover 
(see Fig. 4) when quarks are introduced. The critical temperature is reduced from 
To ~ 270 MeV for the pure gauge system to a transition temperature T c ~ 215 
MeV which almost coincides with the characteristic temperature for the crossover to 
chiral symmetry restoration. In the chiral limit with massless quarks (mo — ► 0), the 
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quark condensate vanishes at T c as a second order phase transition. The crossover 
temperature found in the PNJL model is not far from recent lattice QCD results 
which give T c ~ 200 MeV for two flavors 27 ) and T c = 192 ± 11 MeV for 2 + 1 
flavors 26 ^ . 
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Fig. 3. Chiral (quark) condensate (tptp) and 
Polyakov loop $ as functions of tempera- 
ture at zero quark chemical potential, cal- 
culated using the PNJL model. 14)1 15) The 
inlay shows a comparison of the quark con- 
densate with Nf = 2 lattice data taken 
from ref. 24 ' . 
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Fig. 4. Calculated temperature dependence of 
the Polyakov loop $ in the pure gauge sec- 
tor (dotted curve) using the effective potential 
of Fig. 2, and prediction using the full PNJL 
model (solid curve) 16 ' at zero quark chemical 
potential. Lattice data taken from refs. 25 '' 26 ' 
are shown for comparison. 



4.2. Finite quark chemical potential 

At non-zero quark chemical potential LQCD "data" are deduced from Taylor 
expansions of thermodynamic quantities in powers of j2 = n/T around /i = 0. For 
example, the pressure is expanded as 



p(T,n)=T^Cn(T)fi n , c n (T)- ' i) " (l>/ll) 



n=0 



n\ dfl n 



(4-1) 

ti=0 



We can thus compare the lattice expansion coefficients c n (T) to those calculated in 
the PNJL model. Results for the first few coefficients 02,4,6 are shown in Figs. 5, 6. 
The comparison is evidently quite successful in view of the fact that no fine-tuning of 
parameters has been made, although one must keep in mind that the quark masses 
used in the LQCD simulations are not the small current quark masses used as input 
in the PNJL model. 

It is instructive to examine the quark number density, 



n q (T, M ) = &P ^ ] = 2c 2 yT 2 + 4c 4 ^ + 6c 6 |^ + ... , (4-2) 



and to compare the PNJL result, including Polyakov loop dynamics, with a cal- 
culation using the "classic" NJL model in which quarks propagate freely without 
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Fig. 5. Expansion coefficient cz(T) of Eq.(4-1) 
predicted in the PNJL model 15 - 1 ' 16 ' and com- 
pared to lattice QCD results taken from 
ref. 8 > 



Fig. 6. Same as Fig. 5, for the Taylor expansion 
coefficients ca,(T) and cq(T). The expected 
high-temperture limits (SB) are also indi- 
cated. 



confinement constraints. Fig. 7 clearly demonstrates the important influence of the 
Polyakov loop: it suppresses the quarks as thermodynamically active quasiparti- 
cle degrees of freedom as the transition temperature T c is approached from above. 
Without this restriction, quarks would pile up in the "forbidden" region below T c . 
Further recent results concerning the quark number susceptibility, 

Xq (T, fx) = M) = 2c 2 T 2 + 12 c 4 fJ 2 + 30 ce^ + ... , (4-3) 

have also been analysed 17 ) using the PNJL model which, at this point, can help 
estimating uncertainties induced by truncation of the Taylor expansion in /j,. 




Fig. 7. Temperature dependence of the quark number density (4-2) at fi ~ .1 GeV calculated in 
the PNJL model 14 - 1 (solid curve) and in the NJL model (with <P = 1, dashed curve). Lattice 
QCD results 7 ' using the truncated Taylor expansion (4-2) are shown for comparison. 
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4.3. Phase diagram 

We finally turn to the phase diagram in the (T, y) plane as derived 17 ^ from the 
two-flavor PNJL model. The result is shown in Fig. 8. The solid line between the 
hadronic phase and the quark-gluon phase represents the chiral crossover transition. 
The dashed line marks a first order transition between the phase with spontaneously 
broken chiral symmetry and the (color) superconducting high-density phase in which 
diquarks accumulate to form a non-vanishing (Cooper pair) condensate. The tran- 
sition between this phase and the quark-gluon sector (marked by the dotted line) 
is second order. Further inspection 17 ) reveals that the critical endpoint shown in 
the figure is not a tricritical point; the three transition lines do not have a common 
junction. It turns out in fact that the location of this point is very sensitive to the 
input quark mass and to the presence or absence of Polyakov loop dynamics. For 
example, the temperature at which the gap A disappears is shifted upward by about 
100 MeV when the Polyakov loop is active, as compared to the standard NJL model. 
A more detailed discussion can be found in ref. 17 ) 
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Fig. 8. Phase diagram calculated in the two-flavor PNJL model. 17 ' The solid curve represents the 
chiral crossover region. Dashed curve: first order transition to superconducting high-density 
phase. Dotted curve: second order transition. 



§5. Conclusions and Outlook 

A quasiparticle approach (the PNJL model) encoding the two basic features 
that govern low-energy QCD, chiral symmetry and confinement, has been developed. 
It operates with "order parameter" fields, the chiral condensate and the Polyakov 
loop, coupled through quarks as quasiparticles with dynamically generated masses. 
Despite its extreme simplicity, this approach turns out to be surprisingly successful 
in confrontations with Nf = 2 QCD thermodynamics on the lattice, at least for a 
temperature range up to about twice the critical temperature T c ~ 0.2 GeV. 

Further developments now include extensions beyond mean field theory, 22 ) and 
the generalisation to 2 + 1 flavors with the additional effects of Polyakov loop dy- 
namics, in order to explore the rich variety of superconducting phases. Many more 
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steps are still ahead, such as replacing the notorious NJL cutoff by a running cou- 
pling strength 28 ^ and establishing contacts to the high temperature limit with in- 
corporation of transverse gluons. Last not least, the hadronic matter phase with 
its composite (meson and baryon) degrees of freedom must be approached with an 
improved effective field theory in order to match with a more realistic equation of 
state at high baryon density. 
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